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ABSTRACT 

The time-dependent force acting on a cylinder under- 
going harmonic in-line oscillations in an otherwise steady 
flow was measured for various amplitudes and frequencies of 
oscillation and mean flow velocity. The experiments were 
carried out in a recirculating water tunnel operating as an 
open channel with a free surface at the test section. The 
time-dependent force has been expressed in terms of a mean 
drag coefficient C ^ aand the Fourier-averaged drag and 

inertia coefficients C. and C as functions of the relative 

d m 

amplitude A/D and the frequency paramter D/VT. 
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NOMENCLATURE 



A amplitude of cylinder oscillation 

C d mean drag coefficient 

Fourier-averaged drag coefficient 

C Fourier-averaged inertia coefficient 

m 

D diameter of test cylinder 

F instantaneous total force acting on the test cylinder 
Re Reynolds number (Re = VD/v) ' 

T period of oscillation 

t time 

U instantaneous velocity 

U maximum velocity in a cycle (U = 2 itA/T) 
m m 

V mean velocity of the uniform flow 

e phase angle 

v fluid kinematic viscosity 

p fl ui d dens i ty 
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I. INTRODUCTION 



Elastic structures of one or more degrees of freedom 
can extract energy from the flow about them and can develop 
catastrophic flow-induced vibrations. The understanding of 
this energy-extraction process is of paramount importance 
if one is either to eliminate or minimize it or to design 
the elastic structure such that it can withstand the oscil- 
lations under the contemplated environmental conditions. 

Even a superficial familiarity with the parameters 
involved in this type of phenomenon shows that the investi- 
gation of the flow-induced vibrations or the flow about 
oscillating bodies is more than an extension of the past 
studies on flow about bluff bodies at rest. In fact, the 
complexity of the problem is increased by an order of 
magnitude. In view. of this fact, it is no wonder that the 
past decade has produced either numerous experimental data 
with the goal of obtaining ad-hoc solutions for specific 
problems, while directing little if any attention toward 
elucidation of the underlying mechanics and toward 
generalization of the results for future applications; or 
numerous analytical models most or all of which had nothing 
to do with the motion of the very medium which provided the 
necessary energy to set the body in motion. In fact, the 
fluid mechanics of the flow-induced oscillations became so 
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incidental to the phenomenon that the models dealt essen- 
tially with black-box-induced oscillations. To be sure, 
these models heavily relied on experimental data partly to 
justify their existence and partly to assign numerical values 
to numerous variable constants imbedded in them. The ability 
of these models to scale, i.e., to permit extrapolation, is 
uncertain. Thus, in the final analysis one is not quite 
sure whether one should use the experimental data within the 
range of their application or the curves fitted to them by 
the empirical models. Suffice it to say that the phenomenon 
is far from understood even for the idealized conditions 
encountered in the laboratory without the alarming conse- 
quences of the real ocean environment where fauna and flora, 
ever-changing ocean currents, and temperature gradients add 
further complications to an already complex problem. Surely, 
over design is not the answer but it may, in the next decade 
or so, be the best available tradeoff with failure. 

Mathematical models of flow-induced vibrations of bluff 
bodies and the response of circular cylinders to vortex 
shedding have been aptly described by Parkinson [1] and 
Currie, et al . [2] and will not be repeated here. 

It appears that among the various models considered so 
far, the "wake-oscillator" model of Hartlen and Currie [3] 
attracted more attention among the students of vibration 
analysis. We will not elaborate here on the attempts made 
by others to introduce one or more additional terms into the 
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model originally proposed by Hartlen and Currie for such 
naive attempts produced only more papers and unrealistically 
defined coefficients than sound information toward the 
understanding of the flow-induced oscillations both in-line 
with and transverse to the ambient flow. 

The in-line oscillations has been subjected to very 
little investigation. Chen and Ballengee [4] examined the 
vortex shedding from circular cylinders in an oscillating 
freestream of 3 Hz with A/D form 15 to 1,000 , D/VT = 0.003, 
and the Reynolds numbers up to 40,000, the vortex shedding 
form a circular cylinder responds instantaneously to the 
freestream variations and that "the instantaneous Strouhal 
number stays sensibly constant at 0.2;+ 0.01." This is 
rather expected since the amplitude of oscillations is many 
times that of the cylinder diameter and the flow in the 
absence of significant accelerations exhibits a quasi -steady 
behavi or. 

Hatfield and Morkovin [5] studied the effect of an 
oscillating freestream on the unsteady pressure on a 
circular cylinder (D/VT from 0.15 to 0.25, A/D from 0.05 to 
0.087, and Re = 50,000). They have found that there is no 
significant coupling between the sma 1 1 - amp 1 i tude freestream 
oscillations and the vortex shedding. Their results would 
suggest that the drag coefficient associated with the mean 
flow would essentially remain constant at its steady state 
value. This study, unlike the previous ones, concentrated 
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on the other extreme of the A/D values. In this case, the 
flow fluctuations in the neighborhood of the structure are 
essentially uncorrelated and therefore relatively ineffec- 
tive. Evidently, for values of D/VT from about 0.1 to 0.5 
and A/D from about 0.2 to 1.0 that the interaction of 
oscillations with the body becomes significant and rather 
complex. The study of the two extremes of D/VT and A/D 
does not shed much light on the understanding of oscillating 
flow about flow about bluff bodies. 

Mercier [6] who subjected cylinders to large streamwise 
oscillations found that the averaging drag coefficient 
significantly increases with D/VT and that the rate of 
increase depends on the amplitude to diameter ratio in the 
range of A/D from 0.2 to 3.0, nD/V from 0.1 to infinity, 
and Re from 4,000 to 16,000. 

Davenport [7] subjected bluff bodies (flat, plate, 
circular and triangular cylinders, and lattice truss) to 
small amplitude oscillations in a water flume and evaluated 
the drag and inertia coefficients through the measurement 
of the rate of damping of the amplitude of oscillations of 
the bodies, i.e., without directly measuring the forces 
acting on the bodies. The experiments were carried out in 
the range of D/VT from 0.02 to 2.0. The Reynolds number 
ranged from about 1700 to 8,000. Davenport evaluated C m 
and assuming to be equal to that for the corresponding 

steady flow at the velocity V and ignoring the term 
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2 

involving U . Since C, does not remain constant as evidenced 
m d 

2 

by the present study and U m cannot be neglected for U n ,/V 
larger than about 0.1, Davenport's and C m values are not 
comparable with those presented herein. 

In a related study, Tseng [8] conducted experiments 
with flat plates normal to the stream undergoing freely 
decaying oscillations. His results have shown that the 
presence of the mean flow significantly increases the damping 
force and that the rate of extinction of the oscillation 
increases monotoni cal ly with speed. 

Tanida, Okajima, and Watanabe [9] found that for a 
circular cylinder oscillating parallel to the flow 
(A/D = 0.14, D/VT from zero to 0.5, and for Re = 80 and 
4,000), the vortex synchronization can be observed in a 
range around double the Strouhal frequency, where vortices 
are shed with a frequency half the imposed one. Tanida, 
et alli's results show that the mean drag reaches its maximum 
in the middle of the synchronization range, i.e., D/VT 
between 0.2 and 0.4, and that the sign of is such that no 
energy can be extracted from the fluid to render the oscilla- 
tions unstable. In other words, in-line oscillation is 
stable for the two Reynolds numbers tested. Tanida conjec- 
tured that instability is likely to occur at much higher 
Reynolds numbers. Unfortunately, their values cannot be 
relied upon since the inertial forces were subtracted only 
approximately and, according to the data reported herein. 



are far from being correct. Thus, their values cannot be 
considered sufficiently accurate to assess the stability of 
in-line oscillations. 

Goddard [10] carried out a numerical solution of the 
drag response of a cylinder to streamwise velocity fluctua- 
tions for Re = 40 and nD/V = 0.019, 0.12, and 3.18, and for 
Re = 200 and nD/V = 0.149. He found that for very 1 ow 
frequencies the instantaneous values of drag correspond very 
nearly to the quasi-steady solution and that for higher 
frequencies the drag anticipates the freestream velocity 

maximum. This work cannot be generalized to higher Reynolds 

/ 

numbers since it is based on the Navier-Stokes equations and 
since the diffusion of vorticity in the concentrated vortices 
for Re larger than about 200 is primarily turbulent. 

Numerical or exact solutions of the Navier-Stokes equations 
cannot take into consideration such a turbulent diffusion 
unless modified through the use of an appropriate eddy 
vi scosi ty . 

DECOMPOSITION OF MEASURED FORCES 

The numerical calculations as well as measurements in 
time-dependent flow yield the resultant force as a function 
of time for a given set of numerical values of the independent 
parameters. Thus it is not possible without a suitable 
hypothesis, to express the force both as a function of time 
and remaining independent parameters as one ordinarily would 
in a closed form solution. Such a working hypothesis is 



particularly necessary for the cable strumming problem since 
the results are to be incorporated into the dynamics of the 
cylinder motion in the form of a forcing function. It should 
be stated at the outset that there is, at present, no 
generally accepted hypothesis to decompose the ti me -dependent 
force into suitable components. As it will be seen shortly, 
even the existing hypotheses, such as the so-called Morison's 
equations, are not applicable to periodic flows with a 
non- zero mean velocity. 

Stokes, in a remarkable paper on the motion of pendu- 
lums, showed that the expression for the force on a sphere 
oscillation in an unlimited viscous fluid consists of two 
terms, one involving the acceleration of the sphere and the 
other the velocity. This analysis shows that the inertia 
coefficient is modified because of viscosity and is augmented 
over the theoretical value valid for irrotational flow. The 
drag coefficient associated with the velocity is modified 
because of acceleration, and its value is greater than it 
would be if the sphere were moving with a constant velocity. 

In general , the force experienced by a bluff body at a given 
time depends on the entire history of its acceleration as 
well as the instantaneous velocity and acceleration. Thus, 
the drag coefficient in unsteady flow is not equal to that 
at the same instantaneous velocity in steady flow. Neither 
is the inertia coefficient equal to that found for unseparated 
potential flow. As yet a theoretical analysis of the problem 
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for separated flow is difficult and much of the desired 
information must be obtained both experimentally and 
numerically. In this respect, the experimental studies of 
Morison and his co-workers [11] on the forces on piles due 
to the action of progressive waves have shed considerable 
light on the problem. The forces are divided into two parts, 
one due to the drag, as in the case of flow of constant 
velocity, and the other due to the acceleration or decelera- 
tion of the fluid. This concept necessitates the introduc- 
tion of a drag coefficient C. and an inertia coefficient C 
in the expression for force. In particular if F is the 
force per unit length experienced by a cylinder, then 

F = 0.5 C d p D I U I U + C m pirD 2 /4 dU/dt 

where U and dU/dt represent respectively the undisturbed 
velocity and the acceleration of the fluid. 

On the basis of irrotational flow around the cylinder, 

C m should be equal to 2 (cylinder at rest, the fluid 

accelerating; otherwise C m = 1 ) , and one may suppose that 

the value of should be identical with that applicable to 

a constant velocity. However, numerous experiments show that 

this is not the case and that C. and C show considerable 

d m 

variations from those just cited above. Even though no one 
has suggested a better alternative, the use of the Morison's 
equation gave rise to a great deal of discussion on what 
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values of the two coefficients should be used. Furthermore, 
the importance of the viscosity effect has remained in doubt 
since the experimental evidence published over the said 
period has been quite inconclusive. 

The drag and inertia coefficients obtained from a large 
number of field tests, as compiled by Wiegel [12], show 
extensive scatter whether they are plotted as a function of 
the Reynolds number or the so-called period parameter U^T/D. 

The reasons for the observed scatter of the coefficients C 

m 

and Cj remained largely unknown. The scatter was attributed 
to several reasons or combinations thereof such as the 
irregularity of the ocean waves, free-surface effects, 
inadequacy of the average resistance coefficients to repre- 
sent the actual variation of the nonlinear force, omission 
of some other important parameter which has not been 
incorporated into the analysis, the effect of ocean currents 
on separation, vortex formation, and hence on the forces 
acting on the cylinders, etc. 

The most systematic evaluation of the Fourier-averaged 
drag and inertia coefficients has been made by Keulegan and 
Carpenter [3] through measurements on submerged horizontal 
cylinders and plates in the node of a standing wave, applying 
theoretically derived values for velocities and accelerations. 
Additional measurements have been made by Sarpkaya [14] of 
the in-lines as well as transverse forces acting on cylinders 
and spheres in a sinusoidally oscillating fluid and it was 
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found that the drag coefficient as well as the inertia 
coefficient for a strictly sinusoidally oscillating fluid 
(no mean velocity) is a function of U m T/D and that the 
effect of the Reynolds number is rather secondary and 
certainly obscured by the excellent correlation of the data 
with the period parameter U m T/D. 

On the basis of the above discussion, one would assume 
that Morison's equation would apply equally well to periodic 
flow with a mean velocity where u = V - U cose and that 

and C will have constant, ti me -i nvari ant , Fourier or least- 

m 

squares averages. This, inturn, implies that and C m are 
independent of the associated flow phenomena. There is, 
however, no a priori assurance in the principles of fluid 
mechanics of theory of models that this is, in fact, the 
case. Thus the effect of the combination of a uniform 
current and harmonic oscillations on the time-average and 
oscillatory forces acting on circular cylinders will have 
to be re-examined and the limits of application of the 
Morison's equation be delineated. 

It is a priori evident that both u = -U cose and 
u = V - U m cose yield the same acceleration du/dt. Thus, 
the force in-phase with the acceleration in Morison's 
equation remains unaffected by the presence of the mean flow. 
The results presented herein show that this is not the case. 
Furthermore, the use of the Morison equation as 
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F = 0.5 p C.(tf-U coso) I V-U cose I + C P ttD 2 /4. du/dt 
d m 1 m ' m 



requires that the time-averaged drag force be calculated by 
increasing the force calculated from the steady flow by a 
factor [1 + 0.5(U /7) ]. The results presented herein show 
that such an analysis appreciably underestimates the measured 
mean forces. It suffices to state that the fluid flow 
phenomena for bluff bodies are significantly affected by the 
combination of currents and harmonic oscillations and that 
the results for steady currents alone and oscillations alone 
cannot be combined to yield reliable estimates of forces due 
to both acting together. 

The time-dependent forces in the present study are 
analyzed according to the following three-coefficient 
equation 

F - 0.5C d pDV 2 ♦ C m P^i(-U m cos?it) - 

C d DU m p <T I C0S T Lt jcosyk 
which may be written as, 



F 

0.5pDV 2 



+ C m iT 2 (U m T/D) (D/VT) 2 sin^t - 



C d( U m T / D ^ 2 (°/VT) 2 | cosyk | cos^y-t 
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in which C and C. are given by their Fourier averages as 
m d 



Evidently, , and are functions of VT/D and 

U m T/D or A/D. They may depend also on the Reynolds number 
which does not explicitly appear in the above expression 
because of the assumptions made in the formulation of the 
basic force equation. 

In the foregoing, neither the coefficient is assumed 
to be equal to the steady-state drag coefficient for a 
uniform flow at the constant velocity V, nor C m and are 
assumed to be identical to those obtained for a strictly 
harmonic oscillation. In fact, the results show that = 
(steady) only for U = 0, and and C m are equal to 
those obtained for the harmonic oscillation only for VT/D = 0. 



2 71 

C m = (2U m T/TT 3 D) (F sine) de/( P U^D) 



0 
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II. EXPERIMENTAL APPARATUS AND PROCEDURE 



A. EQUIPMENT 

1 . NPS Water Tunnel 

The experiments were performed in a recirculating 
water tunnel (Figure 1) which had a capacity of approximately 
500 gallons. The galvanized test section was four inches 
wide, eight inches high, and sixteen inches long. A low-rpm, 
high-capacity, fourteen- inch-diameter-d is charge centrifugal 
pump was used to circulate the fluid through the test 
section. The velocity of the fluid was regulated by a 
butterfly valve arrangement which was situated downstream of 
the test specimen. Velocities of .7 to 1.5 fps were obtained 
by adjusting the vains of the butterfly valve. 

2 . Pendulum and Motor Arrangement 

The harmonic motion in this experiment was obtained 
by a small, variable speed, electric motor and crank arrange- 
ment (Figure 2). The crank was attached to the arm of a 
pendulum 41 inches from the pivot point. The overall length 
of the pendulum was 56 inches. The frequency of the pendulum 
oscillation was regulated by the speed of the electric motor. 
The frequency of oscillation varied from 0.6 to 9.0 c.p.s. 

The amplitude of the oscillation was set by adjusting the end 
of the crank on a rotating disk (Figure 2). The amplitude 
"A", ranged from 3/16 to 1.5 inches. 
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3 . Test Specimgn 

Circular cylinders of 3/4 inch and one inch were 
used as the test specimen. The 3/4 inch cylinder was made 
of plexiglass and the one inch cylinder was made of aluminum 
tubing. The test specimen was switched to aluminum tubing 
because the inertial force experienced was less than that 
experienced with the plexiglass cylinder. 

The cylinders were attached to the arm of the 
pendulum by means of an aluminum yoke arrangement (Figure 3). 
The action of the pendulum caused the cylinder to oscillate 
in the center of the test section and in-line with the 
stream flow. 

In order to oscillate the cylinder with the 
pendulum arrangement, it was necessary to operate the tunnel 
with the top of the test section open. The water level was 
maintained as close to the top of the test section as possible 
in order to make the free surface effects negligible. 
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Figure 1. NPS water tunnel. 
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Figure 2 . Electric motor and crank arrangement. 







SCALE 3 ONE HALF DIMENSIONS: INCHES 

Figure 3. Aluminum yoke/force transducer 
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D = 1.0 INCHES 
A/D = 0.5 
T = 0.485 sec 




g 




Figure 4. Force and acceleration traces. 
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D = 1 .0 INCHES 
A/D = 0.5 
T = 0.162 sec 



IN-LINE FORCE 



WET 



DRY 



ACCELERATION 





Figure 5. Force and acceleration traces. 
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Sensor and Sensor Support 



4 . 

The mean velocity of the fluid was derived from 
the pressure differential of a pitot tube located upstream 
of the test specimen. 

The instantaneous acceleration of the pendulum arm 
was monitered by an accelerometer attached to the arm of the 
pendul urn . 

The in-line force- measuring device consisted of a 
cantilever beam as part of the yoke arrangement. Four 
piezo resistive strain gages were mounted on the cantilever 
beam and properly waterproofed. A drawing of this force 
transducer is shown in Figure 3. The transducer was 
calibrated by hanging loads at the mid section of the 
cylinder in the vertical direction. 

B. PROCEDURE 

It was realized, during the early stages of the 
experiment, that a procedure of recording both the wet and 
dry oscillations at a common frequency would be necessary 
for a complete evaluation of the data. The wet oscillating 
force was made up of an inertial force and a resistive force 
produced by the fluid. In order to separate the forces, the 
dry force was recorded at the same frequency as the wet 
force and then subtracted from the wet force to give the pure 
resistive force. 
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To accomplish the physical recording of the wet and 
dry force oscillations, the water level in the tunnel was 
raised and lowered for each frequency tested. The proper 
status of the cylinder was noted on the chart recordings. 

Figures 4 and 5 show the wet and dry forces as recorded 
for two typical frequencies. The acceleration curve was 
used as a reference in matching the wet and dry run curves. 
It was convenient that the acceleration traces were nearly 
perfect sine curves . 

From traces similar to the one shown in Figure 4, the 
in-line force was read and recorded from both the wet and 
dry curves. The result was then punched on IBM computer 
cards for every 1.0 or 2.0 mm's. 

Finally the drag and inertia coefficients were calcula- 
ted through the computer program given in Appendix A. 
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III. DISCUSSION OF RESULTS 



The results of the present investigation will be 
reported in terms of C^, , and C m as a function of D/VT 

or V T/D for various values of A/D. Figures 6 through 9 show 
as a function of D/VT for A/D = 0.25, 0.5, 1.0, and 1.5. 
In all cases is equal to its steady-state value for 
D/VT = 0, i.e, for f = 0. The mean-drag coefficient reaches 
a maximum value at D/VT between 0.2 and 0.35. The value 
of D/VT at which this maximum occurs shifts from 0.35 to 0.2 
as A/D increases. Furthermore, the maximum value of 
increases with increasing A/D. For example, for A/D = 1.5 
and D/VT = 0.1 (a wave of 15 ft amplitude and 20 second 
period superimposed on a current of 5 ft/sec velocity acting 
on a 10 ft diameter cylinder), the mean drag acting on a 
cylinder is increased by nearly 50% 

Upon reaching a maximum, decreases sharply and as 
evidenced by the results of A/D = 0.25, reaches a value 
nearly equal to that corresponding to a steady flow at the 
velocity V. Thus, for a given set of V and D, in-line 
oscillations of either very small or very large frequencies 
do not alter the mean resistance, e.g., for A/D = 0.25, C d 
remains essentially constant in the range 0 . 45 <D/ VT<0 . 1 or 
in the range 2>VT/D>10. In other words, when the frequency 
1 /T = f exceeds twice the Strouhal frequency, there is no 
more vortex synchronization. 
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From a practical point of view, the increase of to 
values several times larger than the steady-state value is 
quite significant. Field studies conducted to determine 
wave forces on cylinders may be strongly affected by the 
presence of mean currents and may yield drag and inertia 
coefficients which are significantly different from those 
obtained with purely harmonic oscillations. It is, of 
course, fully realized that the ocean waves are not usually 
simple harmonic, but are rather irregular. 

Figures 6 through 9 also show two additional facts 
of special importance. Firstly, as far as is concerned, 
thereis no special significance of the case where U m >V. A 
simple calculation shows that U m >V for D/VT>0.637 for 
A/D = 0.25; for D/VT>0.318 for A/D = 0.5; for D/ VT >0 . 1 5 9 
for A/D = 1.0; and D/VT>0.106 for A/D = 1.5. Secondly, 
is not equal to (steady) ( 1 +0 . 5U^/V~^) . In fact, only 
for very small values of U^/V that the time-averaged drag 
force can be calculated by increasing the force calculated 
from the steady flow by a factor ( 1 +0 . 5U^/V”^ ) . 

• The Fourier-averaged drag and inertia coefficients 

and are shown in Figures 10 through 13 and 14 through 17 

respectively as function of VT/D for A/D = 0.25, 0.5, 1.0, 

and 1.5. The coefficients C. and C have also been calculated 

d m 

through the use of the least-squares method and a modified 
least-squares method. Since the results differ very little 
from those obtained through the use of the Fourier analysis. 
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only the Fourier-averaged coefficients are presented herein. 
On each graph, the values corresponding to VT/D = 0 are 
taken from Sarpkaya [14]. 

Firstly, it is noted that both and exhibit larger 

scatter than C.. The reasons for this scatter are rather 
d 

understandable. For a given A, D, and V, A/D is fixed and 
only T is variable in VT/D. Thus, very small values of 
VT/D correspond to very small values of T or to rather high 
frequencies of oscillation. This means that not only the 
fluid forces but also the inertial forces due to the mass of 
the oscillating system are large. Furthermore, at very high 
frequencies, the acceleration is not exactly harmonic and 
the second order accelerations may introduce some error into 
the evaluation of the data. In passing, it should be noted 
that attempts to increase T by decreasing V so as obtain the 
same value of VT/D are not very desirable since this leads 
to smaller Reynolds numbers. This and similar other facts 
simply point out the extreme difficulty of experimentation 
in unsteady flows and are sufficient to put the critical 
reader in touch with his wisdom. 

For large values of VT/D, the frequencies are very low 
and give rise to rather small fluctuating forces about the 
mean drag force. Thus, the decomposition of the small 
alte mating component of the total force into in-phase and 
out-of-phase components, i.e., into drag and inertial forces, 
leads to a number of difficulties. In this range, small 
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errors in phase between, say, zero acceleration and the 
actual position of the cylinder in the cycle lead to 
relatively large errors in the evaluation of and C m . In 
spite of this, however, the data are sufficiently consistent 
to draw quantitative conclusions. It is apparent from 
Figures 10 through 13 that increases rapidly with 
increasing values of VT/D and remains positive throughout 
most of the VT/D values. Only for A/D = 0.25 and 0.50 and 
for very small values of VT/D that becomes negative. 

This range of VT/D values nearly fall in the velocity range 
where U m >V. Evidently, the cylinder cannot extract energy 
from the fluid when (^>0 . Thus, the negative values of 
in a small region of VT/D values imply that it is possible 
to excite in-line oscillations of very small amplitudes at 
frequencies five or more times the Strouhal frequency. 

Such an oscillation has not previously been observed partly 
because no one has excited the elastic test cylinders at 
such high frequencies and partly because the cylinders were 
allowed to oscillate in all directions with the consequence 
that in-line oscillations may have been obscured by rela- 
tively larger transverse oscillations. 

The inertia coefficient C is shown in Figures 14 

m 

through 17 as a function of VT/D for various values of A/D. 
For small values of VT/D, C m approaches its ideal potential- 
flow value of unity. Near Strouhal frequency, however, C m 
decreases sharply, particularly for small amplitudes of 
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oscillation. In other words, near Strouhal frequency, the 
force acting on a cylinder undergoing small amplitude 
in-line oscillations is primarily from drag. This is rather 
interesting if we consider the fact that for small values 
of A/D, say A/D = 0.25 or U m T/D = 1.57, the force acting on 
a cylinder oscillating in a fluid otherwise at rest (i.e., 

V = 0) is essentially of inertial nature [14]. Thus, the 
presence of a mean flow superimposed on harmonic oscilla- 
tions can significantly alter not only the magnitude of 
the inertial force but also the nature of the entire flow 
structure about the body. For large values of VT/D where 
U /V is rather small (U /V = 0.15 for A/D = 0.25, and 
U m /V = 0.30 for A/D = 0.50), the forces acting on the body 
follow the changes in velocity, as noted in connection with 
the discussion of , and the inertial component of the 
force is negligibly small. However, it must be emphasized 
once again that the quasi-steady model is grossly inadequate 
to predict the average drag and the response of the instan- 
taneous force to changes in velocity for large periods or 
for extremely small frequencies bears no relation to the 
existence of an actual quasi-steady state of flow. 
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IV. CONCLUSIONS 



The experimental investigation of the in-line oscilla- 
tions of a circular cylinder in a flow with a mean velocity 
V has yielded the force coefficients , C^, and C^ and has 
shown that: 

(a) The mean flow has significant effects on C^, C^ 
and C^ and that the results of the experiments with harmonic 
oscillations in a fluid otherwise at rest are not applicable 
to oscillations of a cylinder in uniform flow; 

(b) In a region of small values of D/VT where the 
frequency of oscillations is five or more times larger than 
the Strouhal frequency, C ^ is negative and energy may be 
transferred from the fluid to the cylinder. This, in turn, 
may give rise to self-excited in-line oscillations in a 
manner similar to those for the transverse oscillations; and 
that 

( c) The most significant region of the D/VT values in 
which the force coefficient differ significantly from those 
obtained in steady flow or with harmonic oscillations in a 
fluid at rest occur in the vicinity of the Strouhal numbers. 
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Figure 6. Mean drag coefficients versus D/VT for A/D 
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Figure 7. Mean drag coefficients versus D/VT for A/D- 
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Figure 8. Mean drag coefficients versus D/VT for A/D 
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Figure 9. Mean drag coefficients versus D/VT for A/D 
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Figure 10. Fouri er-Averaged drag coefficient versus VT/D for A/D = 0.25. 
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Figure 11. Fouri er-Averaged drag coefficient versus VT/D for A/D 




42 




43 



LO 

CM 

• 

o 

M 

Q 



S® 



0 / 

/ 



/ 



/ 



© / 
/ 



/ 

/© 



/ 

/© 

/<3 



6 



9 

© 

/ 



eg 



Q 

£ 



CT> 



P 



< » 



oo 

o 



o 

CO 

o 



©\ 



V N.° 

<N 



©' 






Nk 

o 



*S. 



@K 



<1 / 

A 



' O 
©!® 

1 <3 



© 



<3 

© ^ 

>-® * 

< 



r. 



® 



< 



o 5^0 

© 

'® 



o 



VO 



CO 



o 



CO 

o 



CO 

« 

o 



•tf- 

o 



CM 

o 



o 



o 

o 



44 



Figure 14. Fourier-Averaged inertia coefficient versus 
VT/D for A/D = 0.25. 
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Figure 15. Fourier-Averaged inertia coefficient versus VT/D 
for A/D = 0.50. 
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Figure 16. Fourier-Averaged inertia coefficient versus VT/D 
for A/D = 1.00. 
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Figure 17. Fourier-Averaged inertia coefficient versus VT/D 
for A/D = 1.50. 
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APPENDIX A 



COMPUTER PROGRAMS 



cccccccccccrcccccccccccccccccccccccccccccccccccccc cccccccccc 

AVERAGE CD AMD CM CALCULATIONS TOP CYLINDERS 
T I f' L = D I * 1 E m 5 ! C N L L S TIME ( T I M FJ / P I CD ) 

BET/=IH 1ENSKJ LLS DI SPLACEMENT ( UMAX* PEP./DI A J 
CL=CYL INNER LENGTH IN FEET 
D I A - D! At ’ E T E P OT CYLINDER IN FEET 
AM P = AMPL I T ODE j) MfTjQfj T*j PEE ‘ r 
P E R - P E R I N D If! SECS ( CHAP T PERIOD/C MARTS PEED) 

Z I = F C P. C E COEFFICIENTS 
CM~ I NE !A C P E F T I C I FNT 
CD = l.i PAG CPErriCITMT 
R E M F = R L M A I N 0 L P. FUNCTION 
A I f a 1 = F r ' IJP I FR C . E FF T c I ENTS 
M= NUMBER n T UAT£ SFT$ 

CN = C COVERT I ON FACTOR 



CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 



D I MENS I ON FORCE (ICO) 

READ IN PAPAMCTEI S WHICH CHARACTERIZE the DATA SET 
N = 12 

DO 100 1=1, N 

READ ( 5 , 1 0 ) 0 I A , AM P, PER,NCARD, UM X , CM , F MM M , NRUN 

R E A D ( 5 , 5 5 ) (FCRCF(K), K= 1 , NC ARD ) 

G = 3 2 . 1 7 A - - 

RHO= 62 . A /G 
PI =3.1 A 159 
CN'U=C. )i 00105 
CL =0.33? 3 334 

V3 AT = ).9 5o.9 

INITIATE DATA SETS 



COMPUTE FORCE Ct EFFICIENTS 

Zl=( 2*Pl:?.«PEP >/(PI*PI*°T*0U*0!A*C L P H 0* A M P ) 
Z2= ( -3~P ER- PER )/ ( 8*r HO*D ! A ^ ° I * C L * A p * A M' 1 ) 

Z 2 L S = ( - t * ° r P P F <• ) / ( ? * P p •- P I * P. H 0 * i ) I A : * C L - A * P * A m P ) 
C COMPUTE R E" A 

BE T .' = 2 * "> I */MP/C! A 

P. E YN 0 = ( Y 3 A p * ■ D I A ) / C N U 

W R IT E ( 6 , 1 5 ) D I A , A f- -ip , P ER , UMX , CM , N RUN 
WRITE (6, 21) 

DO t 5 K = 1 , NCA&D 
SUM =0.0 

SUN = S U 1 + FORCE (K) 

65 corn I NUF ■ 

CORR = SUM/f CARD 
DO 0 5 K=l, NC A RD 
FOPCF(K) = FCPCE(K) - CORR 
FORCE (K) = - FORCE (K) 

85 CONTINUE 

FMM’M = I- MM m - CORR 
• F M m M = - FMMM 

' TIME = 0.0 

CM=0 .0 
C D = 0 . 0 
CMLS---0 
C D L S = 0 
FA A- 0. 0 
FBB = 0 . 0 
FCC= 0.0 
FL'D=0. 0 
FI E = 0 . 0 

DC l.T A T - 1 ,0/f.CARD 
J = NC ARD- 1 
DO 200 K = 1 , J 

F=C.\*( F"f'Cf: ( K ) +l f ROE (K ♦ 1 ) ) /2.C 
ALPHA*.:* PI * ( r I Mi ► Df L T A T / 2 . 0 > 

S INA = S IN ( ALPHA ) 



OOOOOOOOOOOOOOOO 
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C.03A = C' r ( A l PHA ) 

f <: I ; j A = 

fco:a= F^crsA 

F l S = r ■:-Cn? A* ADS ( CP‘. A) 

I LS - F'-'C:?: a * / 3 S ( CO S A ) 

V t A c > - 2 <■' P I*r*l . O^OO^'C r, 3 A 1 ”*" O 3A C'S A’^'C'' SA* r 
rrirtS=2* f> I M «C,. SARAHS (C^SA) *1 .(-U() F 
Ff Co =2* > I *f »S ! -\Y *COSA*AP c (f PSA ) *1 .OCCO*i : 

rnos=2 -p i * c y ' t ! *•:/ ~ s : r a* l.ooo j*f 

r- r. £ C- = 2 *-' I * r »: F * S ! MA«1 .0000*F 

FA A - F A /.j + f A A 

FF 13= FOBS ♦ F p K 

F C C= Ff C j «-f'CC 

F 00 = FUL'S +f urj 

F L E = F E E S ♦ F C F 

CU=FSI V.+CM 

co=i cosa+co 

CMIS=FS!-IA + CMLS 
CDLS=FLS +CI3LS 

V.RITE ( 0, 3 J ) 1 I RE. ALPHA, COSA ,S INA, F , FCGS A, FS IMA 
T IMF =T I 4 E «- D L L T A T 
CCUT jr|U£ 

T ! ME = T I 3F>DL L TA T 

F=c;:*(F.jRr.E ( i )+rn>CE<MCAPD) 1 / 2.0 
Al P IA = 2- R ! y: (TIRL + DE'. T‘ t / 2.0) 

S!»iA=S r; ( ALPHA ) 

CfJ5A = CCS ( ALPHA ) 

F S ! i ; A = F I fJA 

F C C S A = F>-CCSa 

F A A S = 2 ” D ! * r * 1 . 0:' GO* C 0 S A * C ° S A *C 0 S A *C '1 S A *F 

ree$=2*p i *f*f*C‘ sa*/ bs (cos a j *i .r,on *-f 

F C C S = 2 * rj ! * r * c ! f i; • • C J S A* AB S ( SA 1- 1.0 0 00 * F 
FDDS - 2 .-° ! * r *r- 1 M/ * S ! f A* 1 . GOGQ*F 
F E L S = 2 P I * F * f *S! I .A* 1 «0000*F 
F A A= FA A f- + F A A 

FBB = FBtiS + FBB . . 

FCC= FCCE +r C L 
F 0 C= FD Os + F UP 
FEfc-FF ES+FEfc 

cm=fs: rr.+c^ 
co=fc:isa + C P 

CVLS-FSr*’IA + C.MLS 

C DL S = FLS +C Dl S 

KR I 3 E ( 6 » 3 3 ) T I l T , AL PH A , CC SA , S I N A , F , F C CS A , FSINA 
CM = Z l*Cf* -DELTA’ 

CD = /2~Cn *LH LTA’ 

CMLS = 2 1-Cf LS *PFL-TAT 
COLS =Z2LS-CULS U E L T AT 

CM,Ff- = ( Zi/2.0) *(FL E-FAA-FCC*F-3B) / (FDD * F A -FCC-FCC ) 



. 0 * P I ) ) * ( F EE *FCC - F DU*F Be ) / ( F DD *F A t - F ! 



corr = ( -a . o *12 / ( 

WP!’E(6,3f>) 

UR ! T F (0 , AO ) CM , CD , B E T A, PF YfJU r CM LS r CDL S, C?*FF , C OFF 

t NOl L = 0 . 0 • " 

W° IT r ( 6, A5 ) 

TIME =0.0 

00 3 00 K = 1 , f :C A RD 
F = Cf.TFD c CE ( K ) 

THE! A 1 - ( ( 2 . C * P I ) / 360 )*A\'GLE 
C 1 = ( ABS < CCS ( T HL- ■'!))) *C r ' S ( THET A 1 ) 

C 2 = RHL"- ( ( y f ■■ x * * 2 . 01/2.0) !: U r A*Cl 

C3= ( ( P !**.? 1 «0» A* S I f J ( T H E T A 1 ) ) / ( U'1X#PEP.) 

Fl=( CM-C 3-CD*C! ) 

F L S = ( C ML S * C 3 - C Dl S * C 1 ) 

F FCR = ( CM FF ft C3~C OF F rt C 1 ) 

F--F/C2 

R E M F = A B 3 ( F ) - A P S ( F 1 ) 
f MAX = n M‘<*C.M/C2 
P F M F = P F M f / F f 1 ^ X 

P L S = ( A 3 S ( F 1 - A B S ( f L S ) ) / F ' 1 A X — 

PI C : = ( ABS < I ) - AL: ( FI R ) ) / I MAX 

AN Gt Eir,GLlC\D.:(';f:iiRO ir ’ , ' LS,RLS ’ 1 ’ ^• KF0R ' FMAXlCi 
T I M f =T IMF.*- DELTA! 



C-FCC ) 
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3CQ continue 

IOO CO N’T IMUE 

10 FO«MA T ( 3 F3 .A, I8,F8.4 t Fin.6,Fa.4,I7) 

15 F 0 P A T ( 1 L' , 2X,'C!A=' , T 3 . A f ?X , • \ M P= • , F8 . ^ 
1 F 8 . A , 2 X , ' UF X= ' , r C . 4 , 2 X , ' C N = • , ) 8 . A , 2 X , • NT 
20 FORMAT { ’ O’ t?X . ,T ! MF/PCP.' , 7X, ' ALPHA' t 7X , ' 
3* , 2X , ' F ' , 7X, ' r-CL f A • , 5X, ' F S IUA' ) 

30 F V c 0 A T ( '0',j.rl2.A t Fi9.4.2F12.A) 

35 FOR." AT CO* , 7X , « C >'.= 1 , 7X , 1 C 0= • , AX, ’ D!. T A= 1 , 
1 ' C^L S - • » 7X , ' CDl.S= • , 7X, »C OF F= ’ , 7X, • CO FF=* 
AO FQ p M AT ( ' O' , 6F 1 2. A , 2T 12 .4 ) 

A 5 FCR”A T ('O' t 7 X , ' ’ I ML ' , 9 X , • F ' , 9X , • F 1 ' , 9X , • 
1 • RLS « , 12 X, ' FTCF ' , 12X, ' RFFR ' ) 

50 FOR'-' AT ( • O' , 6F12.A, AF15.A ) 

55 FORMAT (F 10. A) 

ST CP 
END 



,?X, • PER*' , 

UN = • ,16) 

COS A • , 15X, 'SIMA 



7 X , • CEY\0= ' , 7 X , 

) 

PEMF • t 9X, ' FL 
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